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Abstract 

The one-electron properties of a certain class of one-dimensional ternary 
quasicrystals are investigated. In particular, we show in detail the presence 
of a special kind of critical states called marginal critical states in these QCs. 
By the use of a real-space renormalization-group method, it is shown that the 
scaling properties of marginal critical states are characterized by stretched 
exponentials. These states are virtually localized, so that their presence may 
make a QC less conductive. 

PACS numbers: 61.44.Br, 64.60.Ak, 71.23.Ft 
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I. INTRODUCTION 



The one-electron properties of the Fibonacci lattice (FL) and the likes have been ex- 
tensively investigated because they are considered to be one- dimensional (ID) models for 
quasicrystals (QCs); these lattices are quasiperiodic and called quasilattices (QLs). One of 
the remarkable features of these QLs is that all the one-electron eigenstates in them are nei- 
ther extended nor localized in the standard meaning, but are of an intermediate kind. More 
precisely, these states exhibit multifractal properties and called critical states. The quan- 
tum dynamics of electrons associated with critical states is characterized by the anomalous 
quantum diffusion, which may be related to the unconventional transport properties of real 
However, it has been revealed recently that most of the ID QLs whose one-electron 
properties have been investigated so far belong to a special class of QLs from the point of 
view of a general classification scheme of QLs, although the class is the simplest and, there- 
fore, the most important.iS Moreover, it has also been reported that there exist new classes 
of QLs whose one-electron properties belong to a different universality class from that of the 
conventional class.! The purpose of this paper is to investigate this subject in more detail. 

Here, we shall summarize briefly the classification of ID QLs.0 ID QLs are in general 
generated by the cut-and-projection methodi from 2D periodic lattices, i.e., the mother lat- 
tices; each of the mother lattices has a hyperscaling symmetry^ characterized by a quadratic 
irrational r which assumes (1 + v5)/2 (the golden ratio), 1 + y/2 (the silver ratio), etc. 
This class of QLs are naturally divided into distinct subclasses, each of which is specified by 
the relevant quadratic irrational. Different QLs in each subclass are specified by different 
windows which are used in the projection method. A QL with a generic window is ternary 
in the sense that there exist three types of lattice spacings, s, m, and I, where the three 
symbols refer to the words, "short", "medium", and "long", respectively. Each subclass is 
further divided into three finer subclasses, type I, type II, and type III, where the sizes of the 
windows for these subclasses are integral, rational, and irrational, respectively, with respect 
to the quadratic field Q[r].Q Of these finer subclasses, the former two are of particular im- 
portance because they have inflation symmetries (self similarities). A generic type I QL for 
each t is represented as a decoration of a binary self-similar QL; that is, all the type I QLs 
for each r form a single mutual- local-derivability class (MLD class) ,0 while the self-similar 
member can be chosen as the representative of the class. The FL and the likes which are 
binary are such representatives. On the other hand, the type II QLs, which are essentially 
ternary, are further classified into an infinite number MLD classes. @ Again, we can choose 
a self-similar QL as the representative of each MLD class, and other members are given as 
appropriate decorations of the representative. 

The quadratic irrational characterizing a type I MLD class including the FL is the golden 
ratio, t = (1 + v5)/2, while those for other type I MLD classes are other "precious ratios 
(means)", e.g., the silver ratio, 1 + Note that there exist an infinite series of type II 
MLD classes which are characterized by a single quadratic irrational. 

Most aperiodic ID structures generated by inflation rules (or, equivalently, substitution 
rules) are not QLs in our definition. The Thue-Morse sequence, the period doubling se- 
quence, and the Rudin-Shapiro sequence are familiar examples. Although our definition of 
QLs is rather restrictive, the one-electron properties of type II QLs remained intact until a 
previous work by the present authors appeared. @ This work will be summarized shortly. It 
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should be emphasized that type II QLs form an overwhelming majority in comparison with 
the conventional QLs, i. e., the type I QLs. 

We can investigate the one-electron properties of a QL by the use of the Schrodinger 
equation in the tight-binding approximation: 

^-1^-1 + ^+1^1 = ^, (1) 

with tjj+i and \L,- being the hopping integral and the probability amplitude, respectively, 
while j runs over all the sites of the relevant QL. We have adopted the so-called bond model 
because we are interested in the properties which are independent of the microscopic details 
of the model. The geometry of a QL can be incorporated into Eq . ([I]) by assuming that 
tjj + i takes one of the three values, t s , t m , and ti, depending on the type of the relevant 
bond. The pioneering work of Kohmoto et al0'0 showed that a similar model based on the 
FL exhibits quite unconventional one-electron properties; that is, 1) the energy spectrum 
is purely singular- continuous or, equivalently, it has zero Lebesgue measure with no point 
component, while 2) all the eigenstates are neither extended nor localized but are critical. 
Both the energy spectrum and the eigenf unctions exhibit multifractal properties. Note that 
critical states are distinctive from localized ones in that the eigenf unctions are not square- 
normalizable. 

The integrated density of states (IDOS) per site, H(E), is a fundamental quantity char- 
acterizing the one-electron properties of a macroscopically homogeneous ID system. The 
scaling exponent a(E) is defined by the asymptotic behavior 

AH = HiE + \e) - H(E - ~e) ~ e a(E \ (2) 



with E being the reference energy and e a small parameter to be extrapolated to zeroJ^M It 
provides the local dimension of the energy spectrum at energy E and satisfies the inequality 
< a < 1. For the purely singular-continuous spectrum of the FL, the exponents a(E) 
distribute continuously within an interval, [a m i n ,a m ax], with < a m \ n < a max < 1. The 
distribution yields the f{a) spectrum of the energy spectruHi.00 The three types of energy 
spectra in one dimension, namely, a point spectrum, an absolutely-continuous one, and a 
singular-continuous one, are associated with localized states, extended states, and critical 
states, respectively. The exponent a(E) gives a measure of the degree of localization of the 
eigenstate at energy level E; a(E) assumes if E is on an isolated energy level, while it 
assumes 1 if E is inside an absolutely-continuous spectrum. Therefore, the spectral char- 
acterization by the scaling exponent a(E) is a basic problem in the studies of one-electron 
properties of QCs. 

In order to analyze the asymptotic behavior of one-electron eigenstates in the FL.several 
renormalization-group (RG) methods have been developed, including the trace map,0'0 the 
exact real-space RGjH and the perturbational real-space RG methods.0 Although all these 
methods are independent in their formulations, they all rely on the inflation symmetry of 
the underlying structure. These methods have also been applied to general ID structures 
with inflation symmetries.&B One of the important consequences of these studies is that 
the one-electron properties of any decoration of the FL are characterized by the common 
kinds of fixed points to those for the FL.0 In other words, all these QLs belong to a single 
universality class with respect to one-electron properties in the standard meaning of the RG 
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theory. It follows then the eigenstates for all the QLs in a single type I MLD class exhibit 
power-law decay with fractional exponents as for the case of the FL. Most of the previous 
generalizations of the FL correspond to this universality class.HH 

Meanwhile, there has also been a large amount of work of mathematicians concerning the 
properties of the energy spectra of ID systems with inflation symmetries.0 The analyses 
have extended over a wide range of ID structures surpassing the list of type I QLs. £3 In 
almost all cases the energy spectra have turned out to be rigorously singular-continuous, 
confirming that the one-electron eigenstates in these self-similar structures are critical. Note, 
however, these results do not specify the detailed properties of the critical states, so that 
there still remains a possibility of different universality classes connected with some self- 
similar structures from that of the type I QLs. 

The present authors have recently reported some preliminary results on the one-electron 
properties of several type II QLs. Indeed, their numerical results on the f(a) spectra suggest 
strongly that all the eigenstates are critical. But they have also found that there exists a 
certain unconventional type of critical states called marginal critical states. With the per- 
turbational real-space RG method, they have proved that the asymptotic behavior of these 
new states is described by a stretched exponential decay, where the corresponding scaling 
exponent a vanishes. The decay is much stronger than power-law decay for conventional 
critical states, and marginal critical states can be thought of as special critical states which 
are almost localized. This result is closely connected with the self-similarity of the type II 
QLs, and it proves that type II QLs belong to different universality classes from that of the 
type I QLs. We present in this paper a full account of the presence of marginal critical states 
in type II QLs based on the exact real-space RG analysis. 

This paper is organized as follows. In Sec.0, we present several examples of type II 
QLs which are discovered by the authors.i Then numerical results for the properties of the 
energy spectrum of a type II QL are presented in Sec. [Ill], and the asymptotic behavior of 
the one-electron eigenstates are analyzed in Sec. [TV] by the use of the exact real-space RG 
scheme. Sec .[V] is devoted to discussions about related subjects and to the conclusion of the 
paper. 



II. TYPE II QLS IN ONE DIMENSION 

We list in Table | representatives of six MLD classes of ID QLs: these MLD classes, 
(A) to (F), have been investigated in our recent work.0 The inflation rule for each of the six 
representative QLs is listed in the last column. Of the six representatives, only the one given 
in the first row is of type I. Note that the first three QLs are obtained from a common mother 
lattice but with different windows. Similarly, the mother lattice used to obtain QL(E) is the 
same as that for the FL, but the window is different. We have used the symbol "QL(E)" to 
denote the QL listed in the fifth row of Table |. Remember that the spacing s is isolated in 
all the QLs in Table | but QL(D). 

In order to illustrate some of the features of these QLs, we shall investigate in more detail 
QL(B). The three spacings of the QL are given by s — 1, m — v^2, and 1 = 1 + a/2 with an 
appropriate unit of length. The inflation rule for this QL is characterized by the inflation 
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whose Frobenius eigenvalue is equal to a = 1 + (= r). The right Frobenius eigenvector 
is proportional to the row vector (s, m, I), while the three components of the left Frobenius 
eigenvector, *(r _1 ,r _2 ,r _1 ), are proportional to the frequencies of the three spacings in the 
QL, where the left superscript "£" stands for the transpose operation. It happens that two 
spacings s and I have an equal frequency; this is the case for QL(C), QL(E), and QL(F) 
as well but not for QL(D). It follows that the average lattice spacing of QL(B) is given by 
a = At- 1 . 

If the inflation operation for QL(B) is repeated n times, we obtain the spacings of the n-th 
generation, s^ n \ m^ n \ and l^ n \ For instance, the first several generations of s^ 1 ' are: = s, 
gC 1 ) = mS; s (2) _ i smSj s ( 3 ) = llsmslsms, • - -. The numbers of the original three types of 
spacings, s, m, and I, included in those of the n-th generation are given by the elements of 
the matrix M n . Using this result, we can show that the total number of spacings included 
in are given by L n = (K n+ i + l)/2, where K n is determined recursively by the equation 
K n+ i = 2K n + K n _i with the initial conditions K — K x — 1. Specifically, L n yields a series 
of numbers: 1, 2, 4, 9, 21, 50, 120, 289, 697, 1682, 4060, • • -. It is evident that L n+1 /L n and 
K n+ i/K n tend both to a = r = 1 + a/2 as n goes to infinity, s^ n \ m^ n \ or provides us 
with the n-th periodic approximants of QL(B). The length of is asymptotically given 
by L n a, so that we need not distinguish it from L n if the unit cell is so chosen that a = 1. 



We investigate the energy spectrum of QL(B). Since there is no analytical method for 
exact evaluation of the energy spectrum, we resort to numerical calculations for periodic 
approximants of the QL. The IDOS obtained numerically is shown in Fig.|l]. Also the 
IDOS of QL(A) is shown for comparison in the same figure. Both the IDOSs exhibit devil's 
staircases because the spectra are purely singular-continuous. 

A singular-continuous energy spectrum includes an infinite number of gaps with varying 
sizes. We shall call an energy level an edge level if it is located on an end of a gap. Energy 
levels accumulate at an edge level from one side only, so that the bottom level and the 
top one are regarded to be edge levels as well. An eigenstate is called an edge state if its 
energy coincides an edge level. Edge states will be found to be key states to understand the 
one-electron properties of type II QLs. 

The multifractal analysis of an energy spectrum can be performed with the formalism 
of Kohmoto.0 The f(a) spectra calculated for approximants of several generations for the 
two kinds of QLs, QL(A) and QL(B), are shown in Fig.[|. The left end of the support of 
the f{oi) spectrum for QL(B) seems to tend, in the thermodynamic limit, to the lower end 
of the a-axis, implying that some of the eigenstates in this QL tend to be localized. This 
is surprising because the distribution of the hopping integrals is narrower in QL(B) than in 



III. NUMERICAL RESULTS 



QL(A). 
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We show in Fig.|] finite-size scaling plots of a m i n and a max . If a quantity obeys power- 
law scaling, the relevant plot exhibits linear behavior. We can obtain the thermodynamic 
limits of these exponents from the plots when they obey power-law scaling. Our estimates 
are: a min = 0.483 and ot max = 0.906 for QL(A). For QL(B), a min does not obey power- 
law scaling and seems to vanish in the thermodynamic limit, while a max exhibits two-cycle 
and power-law scaling seems to be preserved during first several steps, where a max in the 
thermodynamic limit is estimated to be 0.6 ~ 0.7 with only one digit of accuracy. 

Similar numerical analyses for other type II QLs in Table | strongly indicate that the 
vanishing of a m i n is a universal feature of type II QLs, being intrinsically different from the 
case of type I QLs. It is also observed numerically that the energy spectrum of a type II QL 
exhibits non-power-law scaling at an edge level, suggesting that the local scaling exponent 
a for the edge level vanishes. In the next section, the scaling exponent a for an edge level 
of a type II QL is shown to be exactly zero, providing us with a direct proof that a m i n = 0. 



IV. A REAL-SPACE RG ANALYSIS 

The authors have investigated in a previous paper! (to be referred to as FN) the edge 
states of a type II QL by means of a real-space RG method.011 This method is, however, 
based on the perturbation theory, and applies only to the strong modulation regime, where 
\t s \ ^> \t m \ ^ \ti\- In this paper we shall employ an exact real-space RG method which has 
been previously applied to type I QLs. EE This method is based on the Dyson equation 
for the one-electron Green function in real space. Since the method does not rely on any 
approximation, it is applicable to more general situations, including the weak modulation 
regime. 

We shall illustrate the theory with its application to QL(B). The decimation appropriate 
for the bottom energy level of this QL is shown in Fig . f|, where all the sites except those on 
the right side of the type s spacings are decimated; the decimation conforms to the inflation 
rule of the QL. This decimation is applicable to any edge level as well. The relevant one- 
electron Hamiltonian is renormalized at each step of the decimation procedure, and, at the 
n-th step, there appear three types of effective hopping integrals, t$, and £ z , and five 
types of effective site potentials, V^$, Vj^J, V»" , v}" , and V t g , which are associated with 
the five types of local environments of the QL.Q The initial values of these eight-parameters 
are given by the bare parameters, tf> = t x and V$ =0 (x,y = s,m,l). The two sets of 
effective parameters corresponding to two successive steps of the decimation are related to 
each other by 

,/ _ t m t s , tit s , _ t[t m 

ls ~ E-V ms ' m ~ E-VJ l ~ E-X' 

t 2 t 2 
V' = V s i + + ' 



E — V ms E — Vis 



t 2 t 2 
V' = V H - h m 



™. •» ' e - V„ E - V„ 



+ + + (4) 
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V' sl = V 8 i + + 



/2 



E — V ms E — X 

j.2 J-/2 j.2 

V' — V 4- s 4- m 4- m 
ls sm+ E-V ls + E-X + E-V ms 

with X = Vu + $ /(E — Vis). Note that the primed parameter, t' m , participates in the 
right hand sides of some of the equations. This 8D map is transformed with the seven 
dimensionless parameters, 

E — V sm E — V ms E — Vu E — V s i 
o = , c = , a 



t.e ta te t 



S 



E — Vis t m ti 

f = -r> 9 = i~, (5) 



, , „ t t 

into a 7D map: 



a! 



b' 



e — bde + bf 2 g 2 

Tf • 

—b + abe — ef 2 



ef 

, b(c-cde + df 2 g 2 + ef 2 g 2 ) 
f(-ce + f 2 g 2 ) 
_ C e-bcde-f 2 g 2 + bdf 2 g 2 + bef 2 g 2 

f(-ce + f 2 g 2 ) ' W 

i be — abce + cef 2 + abf 2 g 2 — f 4 g 2 



e 



g' 



f(-ce + f 2 g 2 ) 

efg 



-ce + f 2 g 2 ' 



A point of the relevant 7D space is represented by the 7D vector: a = (a, b, c, d, e, /, g). 

The 7D map yields a sensible result only when it has a fixed point. This is the case 
when E is tuned to an edge level. Since the edge levels cannot be estimated analytically, 
we have to obtain them numerically. We take the bottom level as the representative of edge 
levels. Its value can be evaluated accurately enough on the basis of a finite-size scaling 
analysis. The RG flow calculated numerically for the bottom level is shown in Fig. [5] up to 
the 9-th step. We find that the flow converges to a/ = (1,1,1,1,1,0,0), which is a fixed 
point of the 7D map. We can conclude from this result and Eq. (|j) that five quantities of 
the type E — V xy (x, y — s,m, I) are all equal to t s at the fixed point, while t m and U can be 
ignored in comparison with t s . Therefore, the quasiperiodic chains of atoms are separated 
eventually into isolated "atoms" and diatomic "molecules" which are combined through the 
transfer integral t s . That is, the fixed point belongs to the limit of the strong modulation 
regime. This means that RG flow brings the system into the strong modulation regime 
even if the system is, at the beginning, in a weak modulation regime. It follows that the 
perturbational RG theory is applicable in a vicinity of the fixed point, and main conclusions 
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in FN are justified for the weak modulation regime as well. Incidentally, we may conclude 
from Eqs.(||) and (Q) that a ~ d and h ~ e in a close vicinity of the fixed point &j (see 
Fig . [D , since all the denominators in Eq. (0J) are of the same order as t s and both t m and ti 
are negligible in comparison with t s . 

A decimation procedure in the real space RG theory is nothing but a coarse-graining 
in the real space. The same procedure effects as a magnification in the energy space. It 
follows that the energy spectrum has a symmetric triplet structure in a vicinity of each edge 
level; the triplet is composed of three bands. The central band of the triplet is derived 
from isolated "atoms" , while the two satellite bands from isolated "molecules" , which have 
bonding and anti-bonding states. This has been numerically confirmed as shown in Fig.|6| 
(see also Fig. 2 in FN). The relative weights of the three subbands are determined by the 
frequencies of the atoms and the molecules appearing in the system to be r _1 : r~ 2 : r^ 1 in 
the order of them in the spectrum, where r = 1 + \pl. 

Let us consider an asymptotic flow in the vicinity of the fixed point of the present RG. 
In a neighborhood of the fixed point of the 8D map, we can approximately replace by t s the 
denominators of the three equations in the first row of Eq. (01), so that we obtain 



,/ _ t rn t s ,/ _ t(t s ,i _ t\ t s , , 

l s ~ . i l m ~ . i l l ~ .2 ' \ ' 

La La L 



1 s 



That is, the relevant 3D subspace of the 8D space is closed approximately in that neighbor- 
hood. Remember that the structures of the numerators in the right hand sides of the three 
equations above are parallel to those of the inflation rule of QL(B) (cf. Table |). The 3D 
map above is reduced to the 2D map: 




(8) 

whose standard linearization in a neighborhood of the fixed point (0, 0) yields 

(9) 





The 2x2 matrix at the right hand side is nilpotent because its square vanishes, and it 
cannot be diagonalized. This gives rise to a non-scaling behavior of the map at the fixed 
point. If new variables are introduced by / = In / and g = lng, the non-linear map (||) can 
be changed into the linear map: 





M= . (10) 



{',)■ 

The matrix M can now be diagonalized, and the n-th member, (/„, g n ), of the orbit gener- 
ated by the map has the asymptotic expression: /„ « g n -\ oc p n , where p = (1 + v^5)/2 is 
the leading eigenvalue of M; it happens that p coincides the golden ratio. Finally, we obtain 
the asymptotic expressions for / and g: 

fn = g n -i ocexp [-cp n ], (11) 

with c being a positive constant, which depends on the system parameters through the initial 
condition of the original 8D map. 



S 



The ratio of the self similarity of the QL under consideration is equal to a = 1 + y/2, and 
the length scale L = L n characterizing the n-th step of the renormalization transformation 
is proportional to a n . Therefore, the asymptotic expression, Eq.(0), is rewritten as a 
stretched exponential: 



/ n cxexp [-(L/£)1 



(12) 



where £ is a characteristic length depending on the system parameters and the exponent v 
of the stretched exponential is given by 



v = hip/ lner. 



(13) 



To be more specific, we obtain v = ln |(1 + y/E) /In 1 + \/2 rs 0.5460 for QL(B). The 
prediction that the asymptotic flow of the RG is characterized by a stretched exponential is 
confirmed numerically as shown in Fig.^. 

If the linearized form of Eq . (|j) had a diagonalizable transformation matrix, the asymp- 
totic flow would obey power-law, L" 1 /", with a being the scaling exponent for the relevant 
energy level E; this is the case for the fixed point associated with a conventional critical 
state appearing in the case of a type I QL. On the other hand, a disordered system in an An- 
derson localization regime will exhibits an exponential law like exp (— L/£) with £ being the 
localization length of an exponentially localized eigenfunction. The one-electron eigenstates 
characterized by stretched exponentials are called in FN marginal critical states. 



V. DISCUSSIONS AND THE CONCLUSION 

In the preceding section, our analysis has been focused on QL(B) but a similar analysis 
is possible for each of the other four type II QLs in Table | as well. Again, the inflation rule 
provides the proper decimation procedure for the edge states, and the asymptotic behavior 
of the flow is governed by a map reduced into a similar 3D subspace of the 8D parameter 
space. The parallelism between the 3D map and the inflation rule of the relevant QL holds 
as well. The relevant 2D maps and some other quantities of all the type II QLs in Table 
| are listed in Table |I[ The exponents v of the stretched exponentials are not universal 
among different type II QLs. The value of v for QL(D) happens to be the same as that 
for QL(E) and so does that of the scaling parameter a, but this does not mean that the 
two QLs belong to a common universality class with respect to the one-electron properties 
because the 2D maps are distinct. 

It is known that different QLs derived from a common mother lattice but with different 
windows have a common Fourier module which characterizes the positions of the peaks of 
the Dirac measures. It is surprising, therefore, that QL(A), QL(B), and QL(C), having 
a common mother lattice, belong to different universality classes with respect to the one- 
electron properties. Remind in this respect that the mother lattice is common between the 
FL and QL(E) as well. 

The presence of a non-trivial length scale £ makes a marginal critical state lack mul- 
tifractality. This can be observed in the eigenfunction of the bottom state shown in Fig. 
H The whole eigenfunction can be divided into two identical parts, which can be divided 
further into two almost identical parts and so on. This hierarchical structure is connected 



9 



with the structure of the RG near the relevant fixed point; a "molecule" formed by two 
"atoms" of each generation turns out to be an "atom" of the next generation. The weak- 
ening of the interatomic coupling with the generation as a function of the relevant scale in 
real space is given by a stretched exponential, so that the eigenfunction lacks self similarity. 
The coupling virtually vanishes at the generations where the length scale L n satisfies the 
inequality L n ^> £, so that £ gives a measure of the scale of "molecules" which carry virtually 
localized wave functions. The estimated value of £ for the bottom state shown in Fig.|8| is 
£ ~ 3, while the size L of an almost localized molecule in the figure is one for which (L/£)" 
exceeds the order of ten. The presence of marginal critical states makes the left edge, a m i n , 
of the support of the f{a) spectrum extend as far as the lower limit a = 0. Since the whole 
distribution of a is shifted significantly downwards in this case, the one-electron eigenstates 
in a type II QL have strong tendency of localization. This may make a QC less conductive 
if its structure is described by a type II QL. 

We shall investigate the relation between type II QLs and the so called circle sequences 
(CSs)H We begin with introducing CSs. A CS is an infinite string of two digits, and 1. 
It is specified by two fractional parameters, A, ( 6 (0,1), with ( being an irrational, and 
its j-th digit with j G Z is given by — [j( — A] with [■ • •] being the Gaussian bracket. 
The case when ( is a quadratic irrational is important; ( (resp. A) plays a similar role 
for the CS to that of r (resp. the window) for a ID QL, and a class of CSs specified by 
( can be divided into three subclasses, I, II, and III; A of a CS belonging to each class is, 
respectively, integral, rational or irrational with respect to the module Z[£]. A CS has an 
inflation symmetry only when it belongs to subclass I or II. A type I CS is related to a type 
I QL such as the FL. 

It can be shown that a CS with A = 1/2 is symmetric with respect to exchanging the 
two types of digits, so that two digits have an equal frequency in this CS. Moreover, if we 
assume further 1/4 < ( < 1/2, each digit is isolated or paired in the CS but a triplet or 
higher multiplets never appear. These features can be observed, for example, in a CS shown 
in Fig. |, for which A = 1/2 and C = (2 - v / 2)/2. 

Quite recently, the authors have succeeded in revealing a general rule which relates a 
type II CS to a type II QL. To see the relation, we shall tentatively locate the digits of a 
CS consecutively on a periodic ID lattice, and then delete all the lattice points on which 
digit "1" is located. The resulting lattice includes a quasiperiodic array of vacancies. This 
defective lattice is shown under a certain condition to have three types of spacings which 
can be denoted by the symbols s, m, and I according to their lengths. This rule converts a 
CS into a sequence of the three symbols s, m, and I. We can show that the latter sequence 
is isomorphic to the sequence of spacing of a ID QL. A CS can be, conversely, derived from 
a ID QL. We do not, however, present here the proof of these results because a considerable 
space is neccesary. Anyway, four ID QLs, QL(B), QL(C), QL(E), and QL(F), are related 
by this relation to four CSs with A = 1/2 and ( = (2 - V2)/2, V2 - 1, (3 - v / 5)/2, 
and (\/3 — l)/2, respectively, while QL(D) is related to a CS with A = (1 + y/E)/8 and 
C = (a/5 — l)/4. The relation is illustrated for the case of QL(B) in Fig.[|. A proof of the 
general rule relating a type II CS to a type II QL will be presented elsewhere. 

Let CS(B) be the CS which is related to QL(B) by the rule above. Then CS(B) and 
QL(B) will belong to a common universality class with respect to the one-electron properties 
because each of the two quasiperiodic objects is given as a decoration of the other. A tight- 
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binding Hamiltonian of the type of the site model is appropriate for the case of CS(B); 
we may assign two site energies V and V\ to the two digits of the CS. This fact has been 
exploited in FN to investigate one-electron properties of QL(B). We should mention in this 
respect that presence of gaps with unusual scaling properties was for the first time reported 
for CS(E) by LuckH 

In conclusion, we re-emphasize that the most remarkable feature of the one-electron 
properties of type II QLs is the presence of marginal critical states, which belong to a special 
kind of critical states with almost localized nature. The asymptotic behavior of a marginal 
critical state is not characterized by power-law scaling being appropriate for a conventional 
critical state but by a stretched exponential. As a consequence, it can be argued that a type 
II QL is expected to be much less conductive than a type I QL, if one assumes the same 
order of local modulation strength for both QLs; the difference originates from the distinct 
global scaling properties of these two QLs. 
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FIGURES 



FIG. 1. The IDOSs of QL(A) (top) and QL(B) (bottom). The values of the hopping integrals 
taken are indicated in the figures. 
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FIG. 2. The f{a) spectra of the energy spectra for periodic approximants of QL(A) (top) and 
QL(B) (bottom). The hopping integrals as well as the sizes L of the periodic approximants are 
given in the figures, where we mean by a size of a periodic approximant the number of sites in 
its unit cell. f(ot) spectra in the thermodynamic limit are to be obtained from these data with 
the finite-size scaling analysis but we do not show them because the results do not have sufficient 
accuracy. 
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FIG. 3. Finite-size scaling plots of a m j„ and a max for QL(A) and QL(B) with the same values 
of hopping integrals as in Fig. 2. The abscissa is 1/ In L, where L is the number of sites in the unit 
cell. 



FIG. 4. The decimation procedure appropriate for the bottom level of QL(B). 
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FIG. 5. The flow of the parameters a, b, c, d, e, /, g is plotted against n, the generation number, 
for t s = -1.0, t m = -0.8, ti = -0.6. 



FIG. 6. The symmetric triplet structure of the energy spectrum of QL(B), where the hopping 
integrals are taken to be t s = —1.0, t m = —0.5, and ti = —0.3. When the leftmost subband is 
expanded repeatedly, the triplet structure shrinks, since the "atoms" and "molecules" are isolated 
in the relevant limit. 
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FIG. 7. Numerical data of /„ and g n (t s = -1.0, t m = -0.9 ~ -0.998, i, = -0.8 0.996, 

n = 1 ~ 7) are plotted against L n /£ with £ being evaluated through finite-size scaling analyses of 
/ n for fixed system parameters. In the inset, the corresponding data of f n and g n -i are plotted 
against {L n /£) u with ^ = 0.546, confirming that f n = g n -i ~ —{L n /€) u - 



FIG. 8. The eigenfunction of the bottom level of QL(B) in an approximant of L = 1682 with 
= -1.0, t m = -0.8, and t t = -0.6. 
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FIG. 9. Geometrical relation between the circle sequence, with A = 1/2 and £ = (2 - \/2)/2, 
and the relevant type II QL, QL(B). 



18 



TABLES 



TABLE I. A list of ID QLs. The quadratic irrational r characterizing the mother lattice of 
each QL is given in the third column, while the scaling ratio, a, in the fourth column; a always 
takes a power of r. In the last column is given the inflation rule which specifies the self similarities 
of the relevant QL. Identical symbols, s, m, and I, do not necessarily imply the same lengths for 
different QLs. 





type 


T 


ratio a 


inflation rule 


(A) 


I 




T 


s' = m, m' = mms 


(B) 


II 


1 + V2 


T 


s' = ms, m' = Is, I = lis 


(C) 


II 


1 + V2 


T 2 


s' = mms, m' = mmlmms, V = mmlmlmms 


(D) 


II 


2 + ^5 


T 


s' = mss, m' = Islss, I' = Islslss 


(E) 


II 


i(l + V5) 


r 3 


s' = mms, m! = mlms, I' = mlmlms 


(F) 


II 


2 + V3 


T 


s' = ms, m! = mlms, I' = mllms 



TABLE II. The 2D maps and the relevant parameters of the type II QLs listed in Table |[ 



a 2D map p v 

(B) 1+72 /' = <?, 9' = fg W+ 0.5460 

(C) 3 + 2V2 f' = fg, g' = f 2 g 2 + ^3 0.7471 

(D) 2 + ^5 /' = fg, g' = f 2 g 1 + 0.6105 

(E) 2 + ^5 /' = fg 2 , g' = fg 1 + ^2 0.6105 

(F) 2 + ^3 f' = fg,g' = fg ^(3 + ^) 0.7308 
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(A)fc=-1.0, t m =-0.6 




(B)fe=-1.0, t m =-O.S, ft=-0.6 
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